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ABSTRACT 


The  Ris.0  National  Laboratory ,  RosJcilda,  Denmark,  atmos- 
pheric puff  dispersion  model  has  been  tested  for  an  atmos- 
pheric-marine environment.  This  three-dimensional  model 
simulates  the  release  of  Guassian  pollutant  puffs  and  pre- 
dicts their  concentration  as  thay  are  diffused  and  advected 
downwind  by  a  horizontally  homogeneous,  time-dependent  wind. 
Atmospheric  characteristics  such  as  turbulence  intensity, 
potential  temperature  gradient,  buoyant  heat  flux  and  maxi- 
mum mixing  depth  have  been  considered.  Model  predicted  pol- 
lutant concentrations  have  been  compared  to  airborne  sampled 
observations.  The  effect  of  coastal  turbulence  not  observed 
by  the  single  point  meteorological  measuraments  made  onboard 
ship  greatly  affects  the  advection  and  diffusion  of  a  plume 
as    it     moves  onshore.  Additional   measurements/predictions 

particular  to  the  coastal  area  will  have  to  be  incorporated 
into  the  model  for  it  to  accurately  predict  the  onshore 
movement    of   pollutants. 
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I.       INTRODUCTION 

The  downwind  transport  and  distribution  of  atmospheric 
pollutants  from  an  isolated  source  over  land  or  water  has 
become  an  important  environmental  factor  in  today's  society. 
The  need  to  understand  the  distribution  of  smoke,  unpleasant 
or  potentially  harmful  foreign  gases  and  perhaps  radioactive 
debris  from  a  nuclear  powerplant  accident  are  becoming  more 
and  more  essential  for  industral  operations  and  construction 
planning.  The  dispersion  of  such  atmospheric  pollutants  is 
commonly  modeled  by  a  standard  Guassian  plume  model  which 
computes   one-hour   average   plume    characteristics. 

The  Meteorology  Section  of  the  Ris0  National  Laboratory, 
Roskilde,  Denmark,  has  recently  developed  a  puff  model  for 
prediction  and  simulation  of  atmospheric  pollutant  diffu- 
sion. The  model  considers  individual  puffs  of  pollutants 
with  specific  release  rates  that  are  advected  by  a  horizon- 
tally homogeneous  wind  over  a  grid.  The  wind  input  may  be 
either  the  measured  wind  from  a  single  point,  a  spatial 
average  or  a  wind  simulation.  The  model  simulates  the 
instantaneous  plume  characteristics  by  adding  a  group  of 
puffs,  growing  in  size,  as  they  advect  with  the  wind.  A 
Guassian    plume   model,        on   the  other   hand,         provides    a   time 
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averaged  concentration  pattern  based  on  a  single  time  aver- 
age wind  vector.  In  the  puff  model,  the  plume  advects  with 
a  time  series  of  actual  wind  data.  Thus,  the  puff  model  is 
able  to  predict  time  varying  concentration  distributions  in 
actual  changing  wind  conditions,  making  it  an  appropriate 
tool  for  dynamical  computations  of  downwind  dispersions  of 
pollutants. 

A  basic  comparison  of  a  puff  model  simulation  and  a  typ- 
ical plume  is  illustrated  in  Fig.  1.  Looking  from  above, 
the  instantaneous  behavior  of  a  plume  being  advected  from  a 
source  by  the  wind  is  shown.  The  outer  cone-shaped  contours 
represent  the  outer  limit  of  the  plume  boundary  and  are 
identical    in   both  Figs.    1     (A)    and    (B)  . 

Fig.  1(A)  shows  an  instantaneous  depiction  of  an  actual 
plume.  The  long-term  average  plume  concentration  is  shown 
on  the  extreme  right  as  a  smooth  curve  with  a  maximum  on  the 
central  axis.  Also  shown  is  tha  instantaneous  plume  concen- 
tration considered  realistic  but  is  of  such  a  short  time 
scale  that   it   cannot   be   predicted   or    easily    measured. 

The  puff  model  prediction  is  depicted  in  Fig.  1 (B) .  The 
circles  show  the  boundaries  of  individual  puffs  of  pollu- 
tants released   from     the   sourcs.         These   puffs      are   advected 
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Figure  1.   Instantaneous  Behavior  of  a  Typical  Plume  and  a 
Series  of  Puffs  from  a  Puff  Model 

and  diffused  downwind  by  a  frequently  updated  wind.  The 
long  term  average  concentration  prediction  of  the  puff  model 
is  expected  to  be  identical  to  the  long  term  concentration 
of  Fig.  1  (A) .  The  short  term  average  pollutant  prediction, 
a  Gaussian  curve  shown  on  the  extreme  right,  is  not  com- 
pletely realistic  but  is  a  reasonable  approximation  to  the 
instantaneous  plume  concentration  profile. 
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The  purpose  of  this  thesis  is  to  evaluate  and  test  the 
general  characteristics  and  capabilities  of  the  Ris0  puff 
model.  The  adjustable  input  data  will  be  varied  and  pre- 
dicted results  for  various  input  combinations  compared.  A 
preliminary  comparison  will  also  be  made  between  predicted 
results  and  data  collected  from  a  coastal  region  using 
observed  meteorological  forcing  data  as  model  input. 


13 


II,       RISg   PtJFF   MODEL 

A.       GENERAL   CHARACTERISTICS 

The  Ristf  Puff  Model  is  a  three-dimensional  computer 
model  used  for  the  prediction  and/or  simulation  of  the  dif- 
fusion and  advection  of  atmospheric  pollutants.  The  puff 
model  technique  is  to  simulate  a  plume  with  Guassian  shaped 
puffs  with  specified  release  rates  within  a  specified  grid. 
The  initial  size  of  the  puffs  is  normally  one  meter  in  diam- 
eter although  this  can  be  easily  adjusted.  The  amount  of 
material  in  a  puff  is  the  release  rate  times  the  elapsed 
time  between  puffs.  Therefore,  a  long  elapsed  time  between 
puff  releases  results  in  a  higher  initial  puff  pollutant 
concentration  than  a  short  time  interval.  This  should  not 
normally  be  of  concern  if  an  adequate  balance  is  maintained 
between   grid    size,    advection   speed  and   puff   release  rate. 

The  location  of  the  puffs  on  the  grid  is  determined  by 
computing  their  movement  for  a  finite  time  step  using  a 
measured  wind  field.  The  growth  and  buoyancy  of  the  puffs 
are  computed  from  simultaneous  specifications  of  atmospheric 
turbulence      intensity   and      stability   and      from   buoyant      heat 
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flux  at  the  source.  An  inversion  cap  through  which  pollu- 
tants cannot  pass  and  the  source  height  where  pollutants  are 
released  are  variable  and  can  easily  be  adjusted.  Grid  dis- 
tances within  the  model  may  vary  from  meters  to  kilometers 
and  time  durations  from  seconds  to  hours  are  possible. 

This  puff  model  has  the  capability  of  monitoring  a  maxi- 
mum cf  twenty-five  sources  of  puffs  and  its  grid  may  contain 
up  to  100  puffs.  A  puff  source  can  be  located  anywhere  on 
the  grid  and  have  a  unique  release  rate,  start  and  stop  of 
release  time,  and  heat  production.  When  the  center  of  a 
puff  moves  outside  the  boundaries  of  the  grid  (either  hori- 
zontally or  vertically) ,  that  particular  puff  is  dropped 
from  memory.  In  this  way  the  model  does  not  store  irrele- 
vant puff  information,  thus  keeping  computer  memory  require- 
ments to  a  minimum. 

A  variable  to  control  the  amount  of  reflection/absorp- 
tion of  the  pollutant  by  the  surface  is  easily  adjusted  in 
the  puff  model.  Such  a  parameter  is  of  great  value  both  in 
actual  dispersion  problems  and  also  for  gaining  understand- 
ings of  the  plume/ surface  relationship. 

The  model  calculates  the  concentration  at  each  grid 
point  by   summing  the  contributions   from  surrounding  puffs 


15 


for  each  advection  step.  The  grid  concentrations  can 
be  allowed  to  accumulate  or  simply  be  updated  with  the  lat- 
est instantaneous  value.  A  minimum  grid  concentration  of 
interest  can  be  set  to  reduce  computer  run  time  by  dropping 
concentrations  too  small  to  be  of  interest. 

The  output  of  the  model  contains  periodic  results  of 
puff  locations  and  concentrations  as  well  as  initial  input 
verification.  The  time  interval  for  the  periodic  results  is 
adjusted  by  the  input  data.  This  recurrent  lineprinter  out- 
put contains: 

•  X-Y  plane  plots  showing  the  position 
of  the  sources  and  of  puffs  inside  the 
grid, 

•  X-Z  plane  plots  of  puff  positions 
for  evaluating  plume  rise  for  aash 
vertical  level  of  interest,  and 

•  a  table  listing  of  the  grid  point 
concentrations  for  each  level. 

A  computer  drawn  contour  chart  of  the  magnitudes  of  the  pol- 
lutant concentrations  is  also  available. 
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When  considering  distance  between  gridpoints  (delta 
X,Y,Z),  only  spatial  resolution  and  computer  resources  need 
be    considered.  Calculated  concentration      accuracy   is      not 

related  to     the   grid-point    separation.  To    insure     that    no 

essential  information  on  individual  puffs  is  "hidden" 
between  grid  points,  the  grid  separation  should  be  adjusted 
dependent  upon  the  size  of  the  puffs  at  the  downwind  dis- 
tance of  interest.  Other  specific  model  configuration  con- 
siderations are  described  in  the  following  sections.  They 
are  also  discussed  in  more  detail  in  the  model  behavior 
chapter. 

B.        HIND    FIELD 

Once  a  puff  is  released,  it  is  advected  based  upon  wind 
data  measurements  at  a  single  point  only,  normally  the 
release  point.  This  limits  the  validity  of  the  model  to 
situations  where  the  wind  field  and  turbulence  can  be 
assumed  to  be  horizonally  homogeneous  throughout  the  grid. 
It  is  therefore  important  to  insure  that  the  data  obtained 
from  such  a  single  point  measurement  is  representative  of 
the  wind  structure  for  the  whole  area  of  interest. 

The  wind  data  are  normally  obtained  in  the  form  of  a 
horizonal  velocity  time  series.   A  vector  sequence  is  formed 
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by  averaging  over  a  convenient  interval.  These  data  are 
read  into  the  model  after  being  segregatsi  into  turbulence 
classes  as  discussed  in  the  next  section. 

C.   TURBULENCE  INTENSITY  AND  DIFFUSION 

The  growth/diffusion  of  a  puff  depends  upon  the  turbu- 
lence intensity.  To  account  for  this  growth,  the  puff  model 
applies  the  theory  for  relative  diffusion  suggested  by  Smith 
and  Hay  (1961). 

The  turbulence  intensity  is  defined  to  be  the  standard 
deviation  (sigma)  of  the  wind  direction  (in  radians) 
squared.  The  sigma  values  are  collected  for  the  same  short 
time  periods  as  the  wind  speed  measurements  used  to  advect 
the  puffs.  Therefore,  the  intensity  of  the  turbulence  which 
governs  the  relative  diffusion  of  the  puffs,  can  be  adjusted 
along  with  the  the  advecting  wind  speed  after  each  time 
step,  if  conditions  warrent. 

A  very  low  value  of  turbulence  intensity  (as  0.0002) 
represents  a  small  standard  deviation  (0.9),  normally  a  sta- 
ble atmosphere  and  a  weak  puff  dispersion/diffusion.  As  the 
atmosphere  becomes  more  unstable,  the  turbulence  intensity 
increases  along  with  an  increase  of  sigma  values  and  plume 
dispersion/diffusion.     While   these   characteristics   are 
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rasentative   of  turbulence  over  land,      they    can   be  applied   to 
over   water  cases   in  a   broad  sense. 

D.  PLUME    RISE 

In  the  vertical  direction,  puff-risa  can  be  accounted 
for  by  Briggs  (1970)  plume  rise  theory.  In  this  case  buoy- 
ancy is  assumed  to  be  conserved  (adiabatic  motion),  and 
pressure  forces,  molecular  viscosity  and  local  density 
changes  are  considered  small  and  are  neglscted.  The  rate  at 
which  a  puff  rises  as  it  is  advected  dowawind  is  a  function 
of  the  buoyancy  flux,  wind  speed,  puff  distance  traveled  and 
stability  of  the  atmosphere.  Plume  rise  is  considered  sep- 
erately   for  each  individual  puff. 

E.  REFLECTION 

The  interaction  of  the  pollutant  with  the  surface  is 
adjustable  and  can  be  easily  changed  in  the  input  data. 
Total  reflection  or  absorption  Dr  a  fraction  between  the  two 
can  be  used. 

F.  LIMIT    OF    MIXING    DEPTH 

The  effect  of  an  atmospheric  lid  (inversion)  can  be 
applied  in  the  model  to  limit  the  vertical  movement  of  the 
pollutant.    The   model  does   not  permit   the  plume  to  rise 
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above  this  cap.  When  a  maximum  mixing  level  is  in  effect, 
it  acts  to  totally  reflect  the  pollutants  in  the  same  manner 
as  total  reflection  at  the  surface.  This  mixing  cap  would 
act  as  an  inversion  when  the  puff  would  be  expected  to  grow 
much  more  readily  in  the  horizonal  than  in  the  vertical 
direction. 
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III.       DATA    COLLECTION 

An  intensive  field  tracer  study  was  psrformed  during  the 
fall  of  1980  and  winter  of  1981  in  the  Saata  Barbara  Channel 
area  of  the  California  coast.  The  work  was  supported  by  the 
Bureau  of  Land  Management  (BLM)  and  performed  by  the  Envior- 
mental  Physics  Group  of  the  Naval  Postgraduate  School  (NPS) 
and  AeroVironment  Inc  (AV)  ,  Los  Angeles,  ZA.  This  study  was 
designed  to  help  validate  and/or  modify  Gaussian  dispersion 
models  for  coastal  use  and  to  build  a  data  base  for  future 
model  development.  Air  pollution  models  in  current  use  have 
not    been    adequately   validated   for  the  over-water   regime. 

In  the  experiment,  a  tracer  gas  (SF6)  was  released  sev- 
eral miles  offshore  from  the  NPS  Research  Vessel  (R/V)  Aca- 
nia.  Ambient  gas  concentrations  as  low  as  10  parts  per 
trillion  (PPT)  were  determined  by  an  array  of  land  based 
sensors,  from  a  small  boat  and  at  various  levels  by  an  air- 
craft equipped  with  a  continuous  SF6  analyser.  A  chart  of 
the  experimental  area  and  locations  of  the  various  platforms 
is    shown    in  Fig.    2. 
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The  aircraft  flew  through  the  plume  at  various  eleva- 
tions offshore  and  overland.  The  plume  transect  tracks  per- 
tinent to  this  study  were  made  parallel  to  the  coast  approx- 
imately one-half  mile  offshore.  The  airborne  sampling 
consisted  of  instaneous  concentrations  (PPT)  at  selected 
points  at  different  levels  over  a  period  of  six  hours.  The 
observations,  recorded  at  locations  10-70  at  altitudes  of  61 
and  91  m  above  MSL  on  January  29,  1980  ars  shown  in  Fig.  3. 
Average  concentrations  over  the  noted  tine  period  are  shown 
at    the  bottom  of  each   altitude  block. 

The  following  marine  meteorological  parameters  were 
measured  onboard  the  R/V  Acania  while  anchored  approximately 
7. U    km  offshore: 


•relative   wind  speed 

•  wind  speed   fluctuation 
•sea   surface   temperature 

•  sky   cloud  cover 

•  relative    wind   direction 


•air  temperature 
•dew  point 
•ship  roll 
•ship  location 


•  inversion   height    (acoustic  sounder) 
•vertical  temperature   and  humidity   profilies 

(shipboard  radiosonde  launch  every   12    hours). 


23 


TRANSECT 
BEG  END 

17:20:50-17:21:14 


17:33:35-17:34:04 


17:20:50-17:59:42 


DATE:  9/29/80 
ALTITUDE:  300  FEET  MSL 


TIME:  1700-1800  PDT 
1/2  MILE  OFFSHORE 


-10 

288   355   488 

56 

142   521   288   421 

156 

372 

6 

305 

471   1634   372 
471    903   1169   39 

6 

-1 

123 

85 

185    628   219 

82 

20 

-45 

-34 

-21 

54    374   506   174 

-26 

30 


40 
DISTANCE  (GRID  COORDINATES) 


TRANSECT 
BEG  END 

17:25:43-17:26:18 


17:37:00-17:37:41 


17:25:43-17:37:41 


DATE:  9/29/80 
ALTITUDE:  200  FEET  MSL 


TIME:  1700-1800  PDT 
1/2  MILE  OFFSHORE 


1342   4967   4967   1753 

707   595   184 

707   3024  5304   3921   520   558   109 

0    2129   7249   6053 

3848  1419   372   166 

-75   335   5081   7922 

2465   634   634   335 

671   3548   6108  3903 

2277  1007   276    93 

316   1679   5192   5821   260   1511   371    317    167 

30 
DISTANCE  (GRID  COORDINATES) 


40 


Figure  3. 


££?£r^S*. Swerved   Plume  Concentrations    (PPIM    at 
Grid   Coordinates  at   5  1    and    91    i    above   MSL.    29 
September    1980. 
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The  tracer  gas  was  released  at  a  fixed  rate  through  the 
exhaust  of  one  of  the  ships  motor  generator  sets.  The  gen- 
erator was  run  at  a  constant  speed  resulting  in  a  constant 
stack  temperature  and  flow  rate. 

The  above  collected  meteorological  information  enabled 
the  data  grid  to  be  established  and  the  atmospheric  wind 
field,  turbulence  intensity,  stability  and  buoyancy  flux  to 
be  derived  for  inputs  into  the  puff  model.  The  puff  model 
prediction  based  on  this  actual  data  formed  the  basis  for 
the  model  performance  evaluatians  carried  out  here.  From 
this  basis,  different  input  variables  were  adjusted  to  note 
the  effect  on  the  advected  concentrations — both  in  relation 
to  each  other  and  to  the  airborne  measurements. 
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IV.   MODEL  BEHAVIOR 

The  Puff  Model  was  run  on  the  NPS  IBM  S/370  Model  3033 
AP  computer  with  two  goals  in  mind;  (1)  familiarization  with 
model  performance  under  actual  conditions  and  (2)  a  compari- 
son of  model  predictions  with  observed  data.  The  two  goals 
are  interrelated  in  the  sense  that  atmospheric  data  col- 
lected in  the  aforementioned  tracer  study  were  used  to  form 
an  initial  prediction  of  the  pLume  dispersion,  and  varia- 
tions of  that  data  were  used  to  evaluate  the  limits  of  the 
model.  The  data  used  as  input  to  the  model  represented  the 
marine  atmospheric  conditions  as  determined  from  R/V  Acania 
meteorological  data  at  the  time  of  the  experiment. 

Proper  grid  spacing  was  arrived  at  by  considering  puff 
spread,  mean  wind  direction  and  the  geographical  area  of 
interest.  With  the  initial  prediction  in  hand,  data  input 
variables  of  the  model  were  adjusted  and  their  effects 
(changes  in  prediction)  noted.  All  model  predictions  were 
compared  with  the  aircraft  observed  data. 

A  7. 4  X  4.3  km  downwind  area  of  intsrast  was  initially 
gridded  into  an  17  X  10  array.  Distances  between  horizonal 
and  vertical  grid  points  were  approximately  435  m  (Fig.  2). 
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Since  many  of  the  aircraft  observation  times  centered 
around  1730  hours  (all  times  are  Pacific  Daylight  Time), 
model  puff  releases  were  initiated  at  1630.  The  30-minute 
wind  speed  averages  obtained  from  data  taken  onboard  the 
ship  between  1630  and  1730  were  4.7  and  4.8  ms~l.  The  first 
advected  puffs  would  be  expected  to  arrive  at  the  back  edge 
of  the  grid  slightly  before  1700,  and  by  1730  a  steady  con- 
sistent plume  would  be  passing  through  the  area  of  the  air- 
craft track.  (The  model  showed,  in  fact,  puffs  leaving  the 
back  edge  of  the  grid  slightly  before  30  minutes  after  puff 
release)  . 

The  average  wind  direction  was  recorded  on  the  ship 
every  15  seconds.  The  standard  deviations  of  the  wind 
direction  (sigma)  were  computed  as  approximate  one  minute 
averages.  These,  in  turn,  were  averaged  over  30  minutes  to 
correspond  with  the  30-minute  wind  speed  averages.  The 
sigma  values  during  the  time  of  interest  were  1.0  and  0.9 
resulting  in  the  very  small  turbulence  intensity  values  of 
.0003  and  .0002. 

A  delta  Z  value  of  33  m  was  used  to  observe  plume  con- 
centrations at  the  altitudes  of  0,  33,  66  and  99  m  above  the 
surface.  These  levels  were  chosed  for  comparison  with  the 
aircraft  transect  altitudes  of  6 1  and  91  m. 
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Fine  scale  vertical  temperature  and  humidity  plots  were 
drawn  based  upon  radiosonde  soundings  taken  onboard  the 
ship.  The  sounding  taken  at  1735  PDT  (Figure  4)  shows  a 
shallow  unstable  layer  near  the  surface  topped  by  an  inver- 
sion extending  to  near  400  m.  A  80  m  dspth  of  the  mixing 
layer  was  subjectively  established.  The  potential  tempera- 
ture  gradient  computed   by   the   formula 


H  -  If  +  -0098  z  (1) 


was    found   to    be    1.0   deg   K/100    m. 

Basic  data  to  determine  source  strength  and  heat  emis- 
sion from  the  ships  stack  were  taken  from  Schacher,  et  al 
(1981).  As   previously      mentioned,      SF6      gas    was      released 

through  the  ships  motor  generator  exhaust  at  a  constant 
rate.  The  stack  temperature  was  250  deg  F,  the  flow  rate 
was  7.13  $  103  m3s-i.  The  SF6  release  rate  was  47.91  lb 
hr-1.  The  top  of  the  ships  stack — considered  to  be  the 
source  elevation —  was  4m.  The  source  strength  was  con- 
verted to   6.04   gm   s- l    for   input    into    the   aodel. 

Heat    emission    (H)    in   KW  was    determined    by   the   formula 
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H    -   AT    *    J=   *    C      *    Flow  (2) 

RT  p 


where 


delta   T  =  temperature    (stack    -  air)    dag    K 
R  =   dry  air    gas   constant 
=    2.87«M0«   erg/g  deg   K 
p  3    103  *   P(mb)    =  dyne/cm2 
Cp  ■    specific  heat    of   dry   air 

=    . 24   cal/g  deg  K 
Flow    a    16.39   *   7.13   *    103  cm 3/s 

The    heat    emission   was  thus   computed  as   15.07   KW. 

Initially  the  model  grid  was  established  after  noting 
the  area  of  maximum  airborne  sampled  concentrations  (between 
points  24  and  43  of  Fig.  3)  and  the  location  of  the  ship. 
It  became  obvious  during  early  model  runs  that,  with  the 
actual  wind  direction  input,  the  model  predicted  plume  was 
being  advected  south  of  the  grid  towards  point  60  on  the 
aircraft  track.  Obviously,  the  steering  wind,  as  measured 
onboard  ship,  was  not  constant  all  the  way  to  the  shore.  A 
northward  turning  of  the  plume  was  indeed  detected  several 
times  during   the   experiments   by    the   aircraft.      To   compensate 
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for  this  effect,  the  source  of  the  plums  release  was  moved 
in  the  model  three  grid  spaces  (  1305  m)  to  the  north  so  that 
the  maximum  predicted  plume  concentrations  would  pass 
through  the  areas  of  the  maximum  airborne  measured  concen- 
trations. No  corrections  were  made  to  the  model  predicted 
plume  concentrations  because  of  this  adjustment.  However, 
one  could  reason  that  the  predicted  concentration  values 
would  be  higher  in  comparison  with  measured  values  since  the 
coastal  turbulence  and  wind  shift — which  would  tend  to  dif- 
fuse the    plume--were   not   considered. 

Initial  model  runs  with  the  small  turbulence  intensity 
classifications  of  .0003  and  .00  02  failed  to  show  plume  con- 
centrations greater  than  1  *  10- l2  in  tha  grid  at  any  level 
other  than  at  the  source.  Apparently,  tha  grid  spacing  was 
too  large  and  the  narrow  plume  was  advacting  between  the 
grid  points.  In  an  effort  to  locate  the  plume,  a  combina- 
tion of  model  runs  were  performed  varying  the  turbulence 
intensity    and  grid  spacing   as   shown  in  Table    I. 

In  this  table,  a  mixing  level  cap  of  80  m  was  in  effect 
for  tha  model  predictions.  No  concentrations  above  that 
level  were  allowed  in  the  computations-  As  previously  men- 
tioned, the  model  mixing  level  cap  totally  reflects  all 
pollutants   back    downward. 
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With  a  grid  spacing  of  435  m  and  a  turbulence  intensity 
of  .05,  a  plume  concentration  covering  four  grid  spaces  at 
the  east  end  of  the  grid  was  produced.  Predicted  plume  con- 
centrations slowly  decreased  as  the  turbulence  intensity  was 
reduced  to  .02  (atmospheric  stability  increased).  At  an 
intensity  of  .01,  the  concentration  dropped  by  about  five 
orders  of  magnitude.  Normally,  one  would  expect  increased 
concentrations  with  increased  atmospheric  stability.  Per- 
haps the  result  noted  here  is  due  to  the  plume  shrinking 
away  from  a  grid  point  (and  becoming  more  concentrated 
between  the  recorded  grids)  with  the  incrsase  in  stability. 

To  increase  the  grid  resolution,  the  grid  spacing  was 
reduced  by  half  to  217.5  m  and  again  to  108.75  m.  With  each 
reduction  the  grid  was  reduced  by  half  in  the  "Y"  direction 
and  doubled  in  the  "X"  direction  thus  keeping  distances 
between  grid  spaces  equal  in  all  directions.  This  of  course 
greatly  increases  the  computational  requirements.  If  only 
plume  predictions  along  the  back  edge  are  needed  (as  in 
Table  I)  ,  the  downwind  grid  distance  may  be  held  constant  at 
435  m  while  the  horizonal  crosswind  resolution  is  increased. 
In  this  way  many  unnecessary  computations  are  not  made. 
However,   the  increased  horizonal  crosswind  resolution   is 
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computed  over  the  entire  downwind  grid,  which  in  this  case, 
is  not  necessary.  A  more  satisfying  solution  to  this  prob- 
lem is  to  install  the  capability  of  using  a  variable  resolu- 
tion grid  with  the  model  so  that  iownwind  areas  of 
particular  interest  can  be  covered  with  a.  dense  grid  while 
other  areas  of  not  so  much  interest  can  be  sparsely  grided. 

In  order  for  the  advection  of  the  plume  to  remain  on  the 
array  when  increasing  the  horizonal  crosswind  resolution  and 
decreasing  the  area  exposed  on  the  grid,  the  plume  source 
was  adjusted  along  the  western  boundary  of  the  grid.  The 
relationship  of  the  vertical  grid  points  to  changes  of  the 
source  location  is  shown  in  Table  II.  The  plume  source  for 
each  grid  resolution  is  noted  with  an  arrowhead.  Grid 
points  that  are  aligned  vertically  in  the  table  have  identi- 
cal locations  and  should  have  the  same  predicted  plume  con- 
centrations. As  mentioned  earlier,  an  increase  of  grid 
resolution  does  not  affect  the  predicted  concentration. 
Notice  that  for  the  same  grid  points  in  Table  I,  the  pre- 
dicted concentrations  with  a  turbulence  intensity  of  .01 
remain  constant  with  changes  in  grid  spacing — only  the  grid 
resolution  was  changed. 
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TABLE    II 

The   Relationship  of   the   Y   Axis    along   the   Western  Grid   Edge 
to  changes  of  Grid  Distance   between   435    and    108.75    m. 


8 
7        8        9 


DELTA 
Y    (m) 

Grid   Points   Along  the 

435 

4 

5                  6                  7 

217.5 

0 

12        3         4        5        6 

108.75 

0123456789 
I 

From  Table  I  it  is  obvious  that  the  435  m  grid  spacing  is 
too  large  and  that  the  higher  resolution  does  indeed  "see" 
concentrations  that  would  otherwise  be  missed. 

The  problem  of  a  increasingly  narrow  distance  covered  on 
the  grid  as  resolution  is  increased  can  sometimes  be  at 
least  partly  corrected  by  reversing  the  IC  and  T  coordinates 
and  adjusting  or  rotating  the  advecting  wind  direction. 
This  can  easily  be  done  with  the  use  of  the  "TURN"  model 
input  parameter.  This  procedure  sometimes  becomes  necessary 
since  one  of  the  grid  directions  is  limited  by  the  width  of 
the  output  printer  paper  to  less  than  or  equal  10  grid 
units. 

The  model  input  variables,  meteorological  and  source 
values  were  adjusted  to  note  their  effect  on  plume  concen- 
trations. A  deeper  understanding  of  how  the  model  works  and 
how  the  atmosphere  affects  dispersion   can  also  be  gained  by 
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such  adjustments.  A  turbulence  intensity  class  of  .05  was 
used,  except  when  studying  intensity  itself,  because  it  had 
previously  demonstrated  a  good  iownwind  grid  coverage  of  the 
plume. 

In  order  to  note  the  effect  of  the  maximum  mixing  level 
on  the  plume  concentrations,  several  model  predictions  were 
run,  varying  only  the  height  to  which  the  plume  was  allowed 
to  rise.  Exact  grid  point  reproductions  were  not  possible 
since  the  model  only  allows  the  height  of  the  mixing  level 
to  be  an  integer  multiple  of  delta  Z.  The  vertical  grid 
spacing  is  therefore  not  equal.  However,  the  antipicated 
trend  of  increased  concentrations  as  the  mixing  level  is 
lowered  is  evident  from  Table  III. 

The  reflection/absorption  of  the  plume  at  the  surface  is 
controlled  by  the  model  variable  "REFLEZ".  Tests  of  the 
extremes  of  total  absorption  (0.0)  and  total  reflection 
(1.0)  were  performed.  The  results  showed  a  50  percent 
reduction  in  plume  concentrations  at  the  west  end  of  the 
grid  with  total  absorption  compared  to  total  reflection  in 
otherwise  identical  model  runs. 

The  model  has  a  self-imposed  limitation  cf  100  puffs 
from  all  sources  on  the  grid.     The  model  will  terminate  if 
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TABLE   III 

Plume   Concentrations   between    Surface   and    99    m   under 
Different    Maximum    Mixing  Levels. 


Max  Mixing 
Level 

(5) 

(6) 

(7) 

(8) 

None 

.44E-8 
.43E-8 
.40E-8 
.36E-8 

.55E-5 
.54E-5 
.50E-5 
.45E-5 

.20E-5 
.19E-5 
.18E-5 
.16E-5 

.22E-9 
.21E-9 
.20E-9 
.18E-9 

Sfc 
33  m 
66  m 
99  m 

80  m 

.57E-8 
.59E-8 
.58E-8 
0 

.7  2E-5 
.74E-5 
.72E-5 
0 

.26E-5 
.27E-5 
.26E-5 
0 

.28E-9 
.29E-9 
.29E-9 
0 

Sfc 
40  m 
80  m 
>  80  m 

30  m 

.64E-8 
.65E-8 
0 

.81E-5 
.81E-5 
0 

.29E-5 
.29E-5 
0 

.32E-9 
.32E-9 
0 

Sfc 
30  m 
>  30  m 

this  number  is  exceeded.  A  balance  must  be  made  between  the 
rate  at  which  puffs  are  released  from  ths  source  (TAD)  and 
the  time  it  takes  the  puffs  to  be  advectad  across  the  grid. 
A  release  rate  of  one  puff  every  40  seconds  was  predomi- 
nantly used  during   this   study. 

The  turbulence  intensity  variable  was  varied  to  include 
conditions  that  are  more  unstable.  As  atmospheric  instabil- 
ity increases,  the  plume  would  be  expected  to  expand 
whereas,  with  stable  conditions,  the  plume  should  remain 
narrow  and  highly  concentrated. 
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With  the  use  of  NPS  contouring  routines  and  the 
subroutine  "DRAW",  a  visual  comparison  of  the  plume  disposi- 
tion and  concentrations  is  available-  Since  the  computed 
plume  concentration  varies  ovar  many  orders  of  magnitude, 
the  concentration  values  were  converted  to  integer  numbers 
by  multiplying  by  1.  *  10 l3  and  then  talcing  the  logarithm. 
These   logarithms  are   then   smoothed.  Thus,      a   contour    plot 

representing  order  of  magnitude  concentrations  was  be  pro- 
duced. As  with  the  model  variables  "MAPTIM"  and  "KPLANS" 
which  control  the  frequency  and  vertical  levels  of  printer 
plots,  "DRAW"  can  be  called  to  contour  concentrations  at  any 
time   period  and    for  any    level   required. 

Plume  concentration  distributions  foe  turbulence  inten- 
sities of  .05,  .10  and  .25,  all  other  variables  constant, 
are  shown  in  Figs.  5,  6  and  7.  As  expected,  the  plume 
becomes  wider  and  the  concentration  decreases  as  the  turbu- 
lence intensity/diffusion  increases  and  the  atmosphere 
becomes   less   stable. 

In  Figs.  5-7,  the  plume  source  was  located  at  grid  point 
(0,6).  The  wide  plume  in  the  Lower  part  of  the  plot  is  not 
real  but  is  a  function  of  ths  smoothing  routine  spreading 
the    early      puffs   more     than   would     be   expected.  Since  the 
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Figure  5.   Orders  of  Magnitude  of  Plume  Concentration  with 
Turbulence  Intensity  equal  to  .05. 
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Figure  6.   Same  as  Figure  5  except  Turbulance  Intensity 
egual  to  -10. 
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Figure  7.   Same  as  Figure  5  sxcspt  Turbulsnce  Intensity 
egual  to  .25.  J 
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smoothing  routine  would  tend  to  smooth  strong  concentrations 
near  the  source,  the  smoothing  should  ba  eliminated  if  the 
primary  interest  is  near  the  source.  Actually  one  would 
expect  the  puffs  to  behave  as  in  Fig.  3r  from  MikJcelsen 
(1979),  showing  the  relationship  between  the  puff  size  and 
concentration,  the  rate  of  puff  release  (TAG)  and  the 
advecting  wind  speed  0. 


A  PUFF  SIZE 


()    Q    e  C 


ACONCENTRATION 


'T* 


TAU-U 


Figure  8.      Relationship   between    Puff    Size,    Concentration, 

Puff    Release   Rate    (rAO)    and   Advecting   Wind    Speed 
0    (Mikkelsen,1979)  .  * 
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Puffs  would  have  to  travel  the  distance  Xmin  before  they 
expand  to  a  size  where  they  effectively  overlap  and  form  a 
solid  plume.  From  Figs,  5-7,  Dne  can  see  that  plume  concen- 
trations have  increased  with  distance  and  that  a  Xmin  has 
been  reached  in  the  middle  to  upper  part  of  the  grid.  It  is 
at  this  point  that  the  puff  modal  would  be  expected  to  accu- 
rately predict  plume  concentrations.  If  the  area  of  inter- 
est is  before  the  present  Xmin,  the  release  rate  of  puffs 
would  need  to  be  increased  so  that  the  successive  puffs 
would  overlap  sooner.  Also  noted  that  as  the  turbulence 
intensity  increases,  the  area  of  maximum  concentration  of 
the  plume  expands  while  the  central  concentration  decreases. 
This  agrees  with  conservation  of  mass  theorey. 

To  appreciate  the  relative  importance  of  the  source 
strength  and  buoyant  heat  flux,  these  variables  (discussed 
in  Chapter  IV)  were  doubled  separately  and  together  and  the 
concentrations  compared  to  the  concentrations  from  the 
actual  conditions.  Little  or  no  changes  in  concentration 
were  noted  when  rhe  buoyant  heat  flux  was  doubled  and  source 
strength  remained  the  same.  However,  when  the  source 
strength  was  doubled  and  heat  flux  held  constant,  the  grid 
concentrations  doubled  as  expected.    Thus,   under  existing 
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conditions,  the  source  strength  was  critical  to  the  pre- 
dicted plume  concentrations  while  the  buoyant  heat  flux, 
within  the  range  tested,  was  not  relevant.  The  vertical 
printer  plots  did  show  an  initial  puff  rise  soon  after 
release  due  to  the  initial  heat  release  but  as  the  puff  rose 
and  expanded,  it  soon  reached  the  ambient  temperature  and 
leveled  off.  The  buoyant  heat  flux  would  probably  be  more 
important  when  dealing  with  a  smaller  scale  grid  or  greater 
heat  release. 
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V.   DATA  COMPARISON 

No  attempt  was  made  to  compare  actual  puff  model  concen- 
tration predictions  at  exact  grid  points  to  aircraft  obser- 
vations for  the  following  reasons: 

•  The  aircraft  locations  were  approximations — the  exact 
locations  were  not  known.  Large  differences  in  predicted 
concentrations  are  seen  with  small  grid  separations  as 
evidenced  in  Table  II. 

•  As  noted  in  Fig.  2,    the  aircraft  obssrvations  were  taken 
over  a  period  of  time  at  different  lsvels — while  the  puff 
model  produced  multilevel  instantaneous  predictions. 

•  As  mentioned  earlier,  the  actual  wind  was  not  constant 
between  the  ship  observation  site  and  the  oppsite  side 
of  the  grid  near  shore.   Since  the  model  advects  the 
puffs  based  upon  ship  observed  wind,  the  behavior  of 
actual  plume  would  be  different  from  predicted. 

•  Calabration  procedures  for  the  SF6  continuous  analyzer 
mounted  onboard  the  aircraft  were  not  available  for 
instantaneous  concentrations  greater  than  1010  PPT. 
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Therefore,  a  question  of  actual  levels  of  SF6 
concentration  in  the  higher  ranges  exists. 
The  puff  model  predicted  concentrations  are  expressed  in 
g/m3  while  the  aircraft  observations  are  shown  as  the  volume 
of  SF6per  unit  volume  of  air  in  PPT-    h   conversion  between 
the  predicted  and  observed  concentrations  was  obtained  by 
computing  the   partial  pressure  and  molecular  weight   of  SF6 
at  standard  pressure  and  temperature.   &  conversion  of 
g/m3  ■  (.63$10-11)  $  observed  concentration  (PPT) 
was  thus  found. 

Generally,  the  aircraft  sampled  concentrations  (Fig.  2) 
show  values  between  100  and  80  00  PPT.  Converting  these 
observed  concentrations  to  predicted  concentration  units 
gives  values  between  .eS^IO-9  and  .50*10-7  g/m3.  These 
observed  concentrations  are  much  smaller  than  the  values 
shown  in  Table  II.  Perhaps  this  difference  could  be 
explained  by  the  fact  that  the  puff  model  advected  the  plume 
toward  the  coast  in  the  same  direction  under  the  same  very 
stable  conditions  as  observed  on  the  ship.  Any  considera- 
tion of  increased  turbulence  and  wind  shifts  near  shore 
would  be  expected  to  reduce  the  actual  plume  concentrations 
toward  the  observed  concentration  levels. 


46 


Increasing  the  turbulence  intensity  to  0.25  and  keeping 
all  ether  variables  constant,  the  concentration  values  would 
decrease  to  the  order  of  magnitude  of  13-* — closer  to  the 
observed  concentrations,  (This  would  require  the  wind 
direction  standard  deviation  to  increase  from  1  to  28). 
However  the  increased  instability  would  :ause  the  plume  to 
spread  over  a  much  greater  area  (Fig.  7)  than  observed  by 
the  aircraft. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  puff  model  has  been  demonstrated  to  be  a  versatile 
working  dispersion  model.  Different  combinations  of  input 
variables  showed  the  expected  reasonable  results.  The  dif- 
ferences between  model  predicted  and  aircraft  observed  plume 
concentrations  do  not  seem  to  be  the  fault  of  the  model  but 
mainly  that  of  the  highly  variable  meteorological  conditions 
found  along  a  coast. 

Probably  the  most  obvious  conclusion  reached  from  this 
study  is  that  predicting  the  behavior  of  a  plume  moving  over 
a  marine  environment  onto  a  coastal  region  has  significant 
problems.  In  all  probability,  atmospheric  boundary  layer 
conditions  7.4  km  offshore  can  be  very  different  from  those 
observed  in  the  more  turbulent  coastal  region.  The  single 
point  meteorological  measurement  at  the  source  should  not  be 
expected  to  adeguately  represent  plume  characteristics  as  it 
nears  a  meteorologically  variable  coastline.  Additional 
observations  (primarily  wind  speed  and  direction),  or  other 
means  of  predicting  the  coastal  meteorological  conditions, 
would  have  to  be  incorporated  into  the  puff  model  to  ade- 
quately handle  this  problem. 
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The  advantage  of  incorporating  variable  grid  spacing 
within  the  puff  model  and  the  obvious  benifits  have  already 
been  discussed. 

Presently,  the  mixing  cap  of  the  puff  model  is  required 
to  be  located  at  an  integer  multiple  of  dalta  Z.  More  flex- 
iability  in  this  parameter  to  include  any  level,  regardless 
of  delta  Z,  would  be  benifical. 

Along  with  the  puff  locations  shown  on  the  lineprinter 
output,  a  maximum  concentration  level  of  each  puff  would  be 
helpful. 

In  future  experiments,  several  aircraft  tracks  should  be 
made  further  out  from  the  coast  in  an  attempt  to  avoid  the 
turbulent  coastal  region.  Obsarvations  thus  obtained  in  a 
noncoastal  environment  would  help  to  varify  the  model  pre- 
dictions without  the  coastal  influence. 
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APPENDIX  A 
MAJOR  SECTIONS  OF  THE  POFF  33DEL 

The  Ris0  Puff  Model  has  been  described  by  (Mikkelsen, 
1979).  The  code  also  is  well  documented  with  comment  state- 
ments- tfith  that  information  and  the  outline  to  be  provided 
in  this  and  the  following  appendices,  the  computational  and 
input/output  procedures  will  be  obvious. 

The  program  and  input  data  are  stored  on  cards  for  the 
sake  of  permanency.  For  efficient  operational  execution, 
the  program  and  input  data  cards  are  read  on  a  disk  within 
the  computer.  The  model  can  then  be  ran  at  will  without 
reference  to  the  original  data  cards.  Minor  changes  can 
easily  be  made  directly  on  the  disk  both  to  the  model  and/or 
data  before  each  execution. 

The  model  can  be  separated  into  the  following  main  sec- 
tions: 

A.  Input  Data 

B.  Initial 

C.  Calculating 

D.  Output 

E.  Error  Diagnostics 

F.  Subroutines 
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These  will   be  described  separately  in  the  following  sec- 
tions. 

A.   INPUT  DATA  SECTION 

The  input  data  includes  the  variables  shown  in  Table  IV, 

TABLE  IV 
Input  Data  Variables  for  the  Puff  Dispersion  Model, 

Wind  History  Potential  Temperature  Gradient 

Turbulence  Intensity  Buoyant  Heat  Flux 

Grid  Dimensions  Minimum  Concentration  of  Interest 

Mixing  Depth  Reflection  at  Ground  Level 
Source  Locations,  Start/Stop  Time,  Strength,  Heat  Emission 

Number  of  Seconds  between  Advection  Steps 

Number  of  Seconds  between  Printouts/Plots 

Number  of  Seconds  between  Puff  Releases 


The  wind  field  and  stability  class  for  the  current  time  step 
are  read  at  the  start  of  the  calculation  section. 

The  variables  listed  above  are  printed  as  a  input  data 
check  and  a  permanent  record  to  accompany  the  actual  output. 
In  most  cases  the  print  command  can  be  overridden  by  YES/NO 
options. 

B.   INITIAL  SECTION 

Based  upon  the  input  data  from  section  (A) ,  the  initial 
section  specifies  and   initializes  parameters  to  be   used  in 
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the  calculating  section  and  is  passed  only  once  during  exe- 
cution of  the  model.  The  grid  and  some  counters  are  ini- 
tialized. Constants  relating  to  reflectance,  mixing  depth 
and  stability  as  well  as  those  controlling  the  size  of  some 
of  the  loops  within  the  model  are  established.  Parameters 
such  as  number  of  puff  releases  per  second,  number  of  advec- 
tion  steps  per  second  and  number  of  advection  steps  per  puff 
release  are  determined. 

C.  CALCULATION    SECTION 

Using  current  wind  and  stability  class  data  read  at  the 
start  of  the  calculation  section,  the  model  advects  the  puff 
centers  and  calculates  the  growth  rate  and  plume  rise  of  the 
puffs.  It  removes  the  puffs  that  have  lsft  the  grid  (hori- 
zonally  and/or  vertically).  The  predicted  concentration  is 
computed  at  the  grid  points  to  include  pollutants  from  all 
nearby  puffs. 

D.  OUTPUT    SECTION 

For  time  intervals  designated  by  the  input  data,  printer 
plots  of  the  X-Y  and  Y-Z  grid  are  produced.  A  maximum  mix- 
ing   level    is   marked  on   the   Y-Z   grid   if   in   effect. 
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These  plots  include  the  source  location  and  a  trace  of  the 
plume  from  the  release  time  to  the  maptims.  Also  printed  at 
this  interval  is  a  X-Y  table  of  grid  concentrations  for  each 
vertical  level  of  interest.  These  concentrations  can  be 
either  accumulated  or  actual  concentrations  at  the  plot 
time. 

Added  to  the  puff  model  is  a  versatic  plotter  routine  to 
smooth  and  contour  the  grid  magnitude  concentrations  of  the 
above  tables. 

E.  ERROR    DIAGNOSTIC    SECTION 

If  the  model  is  directed  by  the  input  data  beyond  the 
limits  of  the  design  of  the  program,  the  program  is  termi- 
nated by  way  of  the  error  diagnostic  section.  It  prints 
comments  relating  to  commonly  made  input  srrors  enabling  the 
user  to  isolate  problems. 

F.  SUBROUTINES 

The  subroutine  "Sigris"  calculates  tha  puff  size  in  the 
horizonal  and  vertical  directions.  It  also  estimates  plume 
rise  associated  with  pollutent  buoyancy. 

The  subroutines  "Ispace"  and  "Rspace"  are  used  in  the 
framework  of  the  printer  plots. 


53 


The  subroutine  "Draw"  converts  the  plume  concentrations 
to  a  logarithmic  values,  smoothes  and  then  contours  them 
using  NPS  inhouse  contour  subroutines.  The  values  are  con- 
verted to  their  logarithm  values  so  that  the  problem  of  con- 
touring over   many  orders   of  magnitude   is   simplified. 
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APPENDIX  B 
PUFF  MODEL  FLOW  CHART 


INPUT 
SECTION 


READ  -  data  input  variables  (primary,  source, 

intensity/  current  stability) ;  averaging 
timestep  (tav)  for  wind  data 


INITIAL  SECTION   1  Initialize  constants 


READ  &  STORE  all  wind  speed  and  direction 
for  current  intensity  class 


CALCULATION 


SECTION 


next  wind  vector 


RECALL  from  storage  one  wind  vector  for  next  tav 


next  advection  step 


Release  puff,  execute  advection  step, 
calculate  concentration 


<, 


£  *advection 
steps 
.GE.  tav, 


•no 


X 


OUTPUT 
SECTION 


£  advection 
<T  steps   ^*" 
.EQ.  MAPTIM 


.yes. 


yes 


Puff  positions, 
grid  concentra- 
tions, contour 


no 


I  End  of  data  | 
1  Termination! 
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APPENDIX   C 
POFF    MODEL    CONTRACTIONS 


CHEMIN — Minimum   grid  concentration    of   intsrest 

DELX,DELY, DELZ--Distance   in  metars  between    grid   points 

DOSE — Allows   the  concentration   natrix   to   accumulate 

DTDZ — Potential   temperature  gradient    (K/M)     (.GE.    0) 

HEAT — Individual   source    heat   emission    (KW| 

ICOLS — Number  of  columns   in  grid     (.LE-    10) 

INST — Instantaneous  concentration  matrix 

ITIME— Start   time 

JROWS — Number  of  rows   in  grid 

KPLANS — Number  of   vertical    levels  in    grid    (includes   surface) 

MAPTIM — Number  of   seconds    between   printer    plots 

NRELSE — Number  of  seconds   to   stop  of    release 

NRMULT — Number  of  sources    (.LE.    25) 

NTADV — Integer   number  of   seconds    between   advection   steps 
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REFLEC — Reflection   at   ground   level    (0.    none; 1.0    total) 
SOORNR — Number  to   identify    source 
SOORST — Strength  for   individual    source    (gn/s) 
STOPRL — Individual   source   stop  time    (s) 
STRTRL — Individual   source   start    time    (s) 
TAO--Integer   number   of   seconds   between   puff    releases 
TURN — Angle  of  rotation   of   wind    direction 
XSOORCE--X   coordinate  of   source    in  grid   units 
YSOORCE--Y   coordinate  of   source    in  grid   UQits 
ZM--Limited   mixing  depth    (m) 
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